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The developing human brain remains one of the few unsolved mysteries of science. Advancements in 
developmental biology, neuroscience, and medical imaging have brought us closer than ever to understand 
brain development in health and disease. However, the precise role of mechanics throughout this process 
remains underestimated and poorly understood. Here we show that mechanical stretch plays a crucial role in 
brain development. Using the nonlinear field theories of mechanics supplemented by the theory of finite 
growth, we model the human brain as a living system with a morphogenetically growing outer surface and a 
stretch-driven growing inner core. This approach seamlessly integrates the two popular but competing 
hypotheses for cortical folding: axonal tension and differential growth. We calibrate our model using 
magnetic resonance images from very preterm neonates. Our model predicts that deviations in cortical 
growth and thickness induce morphological abnormalities. Using the gyrification index, the ratio between 
the total and exposed surface area, we demonstrate that these abnormalities agree with the classical 
pathologies of lissencephaly and polymicrogyria. Understanding the mechanisms of cortical folding in the 
developing human brain has direct implications in the diagnostics and treatment of neurological disorders, 
including epilepsy, schizophrenia, and autism. 



One of the most intriguing, unanswered questions of the 21st century is what drives the formation of folds 
during human brain development. Fig. 1 illustrates the folding pattern of the human brain, which plays a 
significant role in health and disease: Recent research has, for example, revealed a positive correlation 
between cortical gyrification and intelligence 1 . Beyond that, the surface morphology is associated with malforma- 
tions, like lissencephaly 2 and polymicrogyria 3 , and with neurological disorders, like schizophrenia 4 ' 5 and autism 6 ' 7 . 
These structural abnormalities originate from early in utero development 8 , and may give rise to functional 
pathologies like developmental delay, refractory epilepsy, and cerebral palsy 9 . Rapid advancements in medical 
imaging provide tremendous opportunities to better understand brain development: They enable the quantitative 
classification of malformations 10 and allow us to scan preterm brains in which folding is still incomplete 8 . Yet, 
despite significant progress throughout the past decade, the underlying mechanisms of brain folding remain 
largely unknown. 

Historically, human brain development has been studied using a merely observational approach 11 , mainly on 
the basis of genetics and cytology. Recent evidence suggests that the mechanical microenvironment plays a crucial 
role in controlling pattern formation 12 . In this regulatory process, mechanics and biology are coupled bi-dir- 
ectionally: On the one hand, the biological microstructure changes in response to mechanical loading; On the 
other hand, the mechanical parameters change in response to biological processes 13 . Here we study human brain 
development using a fully coupled biomechanical approach to provide a holistic and predictive understanding of 
the mechanobiology of gyrogenesis. 

The popular hypothesis that axons pull on the cortex to initiate cortical folding 14 has recently been questioned 
and disproven experimentally 15 . Current studies based on magnetic resonance imaging now suggest that differ- 
ential growth 16 - with a faster growing surface layer and slower growing inner core - drives the formation of 
cortical folds 17 . From other organ systems like the airway 18 , the intestine 19 , and the gut 20 , but also from plants 21 and 
seashells 22 , we know that differential growth can generate sufficient compressive residual stresses to induce 
structural instabilities such as buckles, folds, or creases 23 , even in the absence of external loading 24 . It seems 
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Figure 1 | Schematic illustration of the human brain. The cross section of 
the brain illustrates its characteristic convoluted surface with ridges called 
gyri and grooves called sulci. The outer layer of gray matter, which 
forms the cortex, is made up of neuronal cell bodies including the cell 
nuclei. The inner core of white matter, which forms the subcortex, consists 
primarily of neuronal axons. Cortical folding maximizes the surface-to- 
volume ratio to increase the number of nuclei and decrease the relative 
distance between them. 

intuitive to hypothesize that differential growth - with a faster grow- 
ing cortex and slower growing subcortex - also plays a major role 
during brain development. 

In this study, we establish a mechanical model for differential 
growth to simulate cortical folding in the developing human brain. 
The key challenge is to identify evolution equations of growth that 
realistically mimic the underlying biochemical processes 25 ' 26 . In the 
simplest case, growth is purely morphogenetic, evolving in time irre- 
spective of the mechanical environment 11 . In more complex cases, 
growth is mechanically driven, for instance by hypertension in arter- 
ies 27 or by hyperstretch in skin 28 . In our model, the cortex, the outer 
surface of the brain, grows morphogenetically at a constant rate, 
while the subcortex, the inner core, grows in response to overstretch. 
Subcortical growth thus mimics the phenomenon of towed 
growth 29 ' 30 , the chronic elongation of axons when stretched beyond 
their physiological limit 31 ' 32 . 

We hypothesize that morphological abnormalities in the human 
brain originate from misbalanced cortical and subcortical growth. 
Our rationale is that mechanically- driven growth has recently been 
shown to explain gyral wavelengths 33 and stress patterns 15 in the 
developing ferret brain. Based on these findings, we reveal mechan- 
isms that cause abnormal brain folding and correlate our results with 
characteristic pathologic observations in lissencephalic, polymicro- 
gyric, and schizophrenic brains 8 . Identifying which parameters affect 
gyrogenesis and surface morphology can provide a tremendous 
insight into the human brain - how it forms and how it functions. 

Results 

Development of brain surface morphology. In the human brain, 
cortical folding takes place primarily in utero, during 24 and 32 
weeks of gestation. In preterm infants, cortical folding is 
incomplete at birth, which allows us to explore the evolution of 
folding patterns using standard medical imaging techniques. Fig. 2 
contrasts magnetic resonance images and simulated brain slices of 
the developing human brain in very preterm neonates at 27, 29, and 
32 weeks of gestation. The magnetic resonance images (Fig. 2a) were 
acquired as part of a four-year prospective neuroimaging study of 
very preterm neonates after receiving informed, written consent 
from the parents. The study was approved by the Research Ethics 
Committee 8 . The magnetic resonance images reveal the spatio- 
temporal evolution of the characteristic convoluted surface morpho- 
logy. Close ups of the blue rectangular cross section (Fig. 2b) 
illustrate progressive cortical growth associated with an increase in 
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Figure 2 | Development of brain surface morphology, (a) Magnetic 
resonance images of the developing human brain in very preterm neonates 
at 27, 29, and 32 weeks of gestation display the evolution of the 
characteristic convoluted surface morphology. The blue rectangle frames 
the two-dimensional rectangular domain for the simulation, (b) Close ups 
of the simulated domain display a progressive increase in sulcal depth until 
the gyri form first self-contact, (c) Simulated morphologies of developing 
human brain section at 27, 29, and 32 weeks of gestation display a gradual 
increase in cortical volume associated with progressive gyrification. 

sulcal depth. The simulated brain slices (Fig. 2c) are based on a finite 
element model of cortical and subcortical growth. The underlying 
continuum model, its finite element discretization, and the chosen 
parameters are specified in detail in the Methods section. The 
simulated brain slices display the progression of growth: At week 
27, the cortex has increased 65% in volume and shallow sulci begin 
to form; At week 29, the cortex has increased 76% in volume and sulci 
become more prominent; At week 32, the cortex has doubled its 
volume, sulci increase significantly in depth, and gyri begin to 
form first self- contact. The simulations predict the characteristic 
features of the developing human brain: gyrogenesis and 
progressive cortical folding between 27 and 32 weeks of gestation 
(Fig. 2c). 

Effect of cortical growth rate on surface morphology. Cortical 
growth depends on various biological processes like cellular differen- 
tiation, growth, proliferation, and migration. Not surprisingly, it 
varies significantly both in space and time. Fig. 3 illustrates the 
effect of varying cortical growth rates on growth, elastic volume 
stretch, gyral wavelength, and gyrification index as predicted by 
our model. For cortical growth rates below a critical value, subcor- 
tical growth is fast enough to compensate for cortical growth: 
Gyrification is suppressed and the surface remains smooth. Increas- 
ing the cortical growth rate beyond this critical value decreases the 
relative amount of subcortical growth (Fig. 3a), which increases the 
elastic volume stretch (Fig. 3b). Although the wavelength remains 
almost unaffected by this change (Fig. 3c), the reduced relative 
subcortical growth reduces the depth of the individual sulci, which 
decreases the gyrification index (Fig. 3d). In summary, with 
increasing cortical growth rate, folds maintain their length, but 
become more shallow. 

Effect of initial cortical thickness on surface morphology. The 

initial thickness of the cortical surface varies significantly throug- 
hout the human brain. Fig. 4 illustrates the effect of thickness 
variation on growth, elastic volume stretch, gyral wavelength, and 
gyrification index as predicted by our model. Increasing the initial 
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Figure 3 | Effect of cortical growth rate on surface morphology. 

Sensitivity of (a) growth <9, (b) elastic volume stretch / e , (c) gyral 
wavelength 2, and (d) gyrification index GI for varying cortical growth 
rates G c . For cortical growth rates below a critical value, no gyrification 
occurs. Increasing the cortical growth rate beyond this value does not affect 
the gyral wavelength but decreases the gyrification index: Folds maintain 
their length, but become more shallow. 
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Figure 4 | Effect of initial cortical thickness on surface morphology. 

Sensitivity of (a) growth #, (b) elastic volume stretch / e , (c) gyral 
wavelength 2, and (d) gyrification index GI for varying initial cortical 
thickness t c . Increasing the initial cortical thickness increases the gyral 
wavelength and decreases the gyrification index: Folds become longer and 
more shallow. 
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Figure 5 | Pathology and simulation of lissencephaly. (a) Magnetic 
resonance images of a lissencephalic brain display an increase in cortical 
thickness, a decrease in depth and number of cortical gyri and sulci, and a 
decrease in subcortical volume, (b), (c) Simulated morphologies for a 
moderately and severely increased cortical thickness display the formation 
of a single shallow gyrus and the complete suppression of gyrification. 



cortical thickness decreases the amount of growth in both cortex and 
subcortex (Fig. 4a). Reduced cortical growth results in less localized 
deformation patterns, reduced axonal tension in the subcortex, and 
reduced subcortical growth. Increasing the initial cortical thickness 
affects the elastic volume stretch only marginally (Fig. 4b). However, 
it markedly stiffens the cortical surface, which increases the gyral 
wavelength (Fig. 4c), and decreases the gyrification index (Fig. 4d). 
In summary, with increasing initial cortical thickness, folds become 
longer and more shallow. 

Lissencephaly. Lissencephaly is a rare malformation of the human 
brain associated with migration disorders and a poorly convoluted 
cortex 9 . During neuronal migration, a significant fraction of neurons 
fails to reach the outer cortex and remains underneath the cortical 
plate. Neuronal misplacement results in a severe defect of cortical 
connectivity and a markedly thickened cortex. Fig. 5 illustrates the 
pathology and simulation of a lissencephalic brain. Magnetic 
resonance images reveal an increased cortical thickness associated 
with a decrease in depth and number of gyri and sulci, and a 
decreased subcortical volume (Fig. 5a). Simulations of moderately 
and severely increased cortical thicknesses display the characteristic 
morphological abnormalities of lissencephalic brains: the formation 
of a few shallow gyri (Fig. 5b) and the complete suppression of 
gyrification (Fig. 5c). 

Polymicrogyria. Polymicrogyria is a malformation of the human 
brain associated with organization disorders and an overly 
convoluted cortex 9 . At inspection, the polymicrogyric brain surface 
appears smooth with only a few aberrant sulci. Cutting the brain 
exposes the characteristic convoluted structure with an increased 
number of secondary folds. Fig. 6 illustrates the pathology and 
simulation of a polymicrogyric brain. Magnetic resonance images 
reveal a reduced cortical thickness associated with an increase in 
number and a decrease in size of gyri and sulci. (Fig. 6a). 
Simulations of a reduced cortical thickness combined with 
physiological and severely retarded subcortical growth display the 
characteristic morphological abnormalities of mircogyric brains: the 
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Figure 6 | Pathology and simulation of polymicrogyria. (a) Magnetic 
resonance images of a polymicrogyric brain display a reduced cortical 
thickness, an increase in number and a decrease in size of gyri and sulci, (b), 
(c) Simulated morphologies for severely reduced cortical thickness at 
physiological and severely reduced subcortical growth rates display 
increased gyrification without and with formation of secondary folds. 

formation of numerous small gyri and sulci (Fig. 6b) and the 
formation of secondary folds (Fig. 6c). 

Schizophrenia. Schizophrenia is a mental disorder with symptoms 
like delusions, social withdrawal, and disorganized thinking asso- 
ciated with regional structural abnormalities in the brain. Fig. 7 
contrasts the pathology and simulation of a healthy and a 
schizophrenic brain. Surface reconstructions reveal regional 
structural differences in the color-coded regions 34 . In comparison 
with the healthy brain (Fig. 7a), the schizophrenic brain displays a 
locally reduced cortical thickness and a locally increased gyrification 




Figure 7 | Pathology and simulation of schizophrenia, (a), (b) Surface 
reconstructions of healthy and schizophrenic brain display a reduced 
cortical thickness and an increase in number and depth of gyri and sulci in 
the color-coded regions of schizophrenic brains, adopted with permission 
from 34 , (c), (d) Simulated morphologies in blue rectangular cross-section 
through the superior temporal gyrus at physiological and reduced cortical 
thickness display an increase in number and depth of gyri and sulci 
associated with an increased gyrification index. 



index (Fig. 7b). Simulations of a reduced cortical thickness in the blue 
rectangular cross-section of the superior temporal gyrus display the 
characteristic morphological abnormalities of schizophrenic brains: 
the regional increase of the gyrification index from 1 .47 in the healthy 
brain (Fig. 7c) to 1.71 in the schizophrenic brain (Fig. 7d). 

Discussion 

Mechanics plays a critical role in brain development. Brain deve- 
lopment was long considered an exclusively morphogenetic process, 
independent of forces, stress, stretch, or strain. There is now 
increasing evidence that mechanics plays a significantly role in 
regulating brain morphology. Here we have explored the role of 
mechanical stretch in human brain development using the nonli- 
near field theories of mechanics supplemented by the theory of 
finite growth. This allowed us to bridge the scales and correlate 
subcellular and cellular events like axon elongation to clinically 
relevant characteristics like surface area and gyrification indices. 

The competing hypotheses of axonal tension and differential 
growth are not mutually exclusive. Mechanical factors are reco- 
gnized as potential driving force for cortical folding, but contro- 
versies around two competing hypotheses, axonal tension and 
differential growth, have slowed down further progress 35 . Axonal 
tension, a mechanism to bring functionally related units topographi- 
cally close together 14 , disagrees with dissection experiments, which 
confirm the general existence of axonal tension, however, not along 
the predicted directions 15 . Differential growth, a mechanism to 
release residual stresses by surface buckling 16 , agrees with stress 
distributions in dissection experiments, however, not with 
physiologically realistic stiffness ratios 35 . Here we combine both 
mechanisms in a unified model with a morphogenetically growing 
outer surface 36 and stretch- driven growing inner core 28 . This allows 
us to predict both axonal tension, in the form of the elastic volume 
stretch / e , and differential growth, in the form of cortical and 
subcortical volume growth S c and $ s . 

Cortical growth induces extreme subcortical deformation and 
axonal stretch. There is a general agreement that cortical growth is 
almost entirely morphogenetic, independent of mechanical factors 17 . 
We know, however, that rapidly growing biological membranes like 
the cortex may induce mechanical instabilities 24 , associated with 
large deformations of the overall system 36 . In the ferret, cortical 
growth is as rapid as 12.7 mmVday per hemisphere in the third 
week after birth, and even increases to 36.7 mmVday per hemis- 
phere in the fourth week 37 . In humans, cortical growth is on the 
order of $ c = 1.8 in first year after birth, and decays towards 
S c = 1.2 in the second year 38 . Not surprisingly, this creates extreme 
subcortical deformation, which our model captures in the form of the 
elastic volume stretch f. Fig. 3 illustrates this effect through an 
increasing elastic volume stretch f in response to increasing 
cortical growth rates G c . 

Subcortical growth is a result of stretch-driven axon elongation. 

On the cellular level, chronic axonal overstretch activates mechano- 
transduction pathways, which collectively result in a gradual increase 
in axonal length 31 . In the blue whale, axonal growth can be as rapid as 
3 cm/day, in the giraffe up to 2 cm/day. In humans, maximum 
axonal growth in utero can be up to 11 cm/month during peak 
growth around the fifth months of gestation 30 . In vitro, under 
extreme conditions, axons tolerate elongation rates of up to 8 mm/ 
day towards a total growth of 10 cm 39 . On the subcellular level, 
axonal growth is associated with an increase in the number of 
microtubules and neurofilaments, an increase in cytoplasm, and an 
increase in plasma membrane, while the ultrastructural appearance 
and the diameter of the growing axon remain virtually unchanged 39 . 
Axonal growth is a self- regulating mechanism: As the axon grows, it 
reduces its stretch, and thereby regulates its length 30 . These 
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phenomena are collectively known as towed growth, the chronic 
axon elongation in response to mechanical stimuli 29 . While initial 
axon elongation experiments use tension as mechanical stimulus 30 , 
recent experiments prefer to use stretch instead 40 . Since stretch- 
driven growth is purely kinematic, it is easier to control and 
manipulate experimentally 39 . In our model, we suggest stretch 
rather than tension 33 as the driving mechanism for subcortical 
growth. Our model inherently captures the self- regulating 
mechanism of growth by activating axon elongation only if the 
elastic volume stretch f chronically exceeds the physiological limit 

JO 41 

Misbalanced cortical and subcortical growth creates morpholo- 
gical abnormalities. Our results confirm our hypothesis that 
morphological abnormalities in the human brain originate from 
misbalanced cortical and subcortical growth: Abnormally slowly 
growing cortices allow axons to almost instantaneously respond to 
growth-induced subcortical deformation so that no folds emerge 
(Fig. 3, left); Physiologically growing cortices create sufficient 
subcortical deformation to trigger axon elongation and subcortical 
growth (Fig. 3, middle); Abnormally fast growing cortices reduce the 
sulcal depth, decrease the gyrification index (Fig. 3, right), and 
provoke the formation of secondary folds 12 (Fig. 6, right). These 
findings agree with qualitative observations in diseased human 
brains 16 : Lissencephalic brains, characterized by a poorly 
convoluted morphology (Fig. 5a) display reduced growth in the 
cortical layer (Figs. 5b,c); Polymicrogyric brains, characterized by 
an overly convoluted morphology (Fig. 6a) display reduced growth 
in the subcortical layer (Figs. 6b,c). 

Our results further demonstrate that morphological abnormalities 
in the human brain could also originate from misbalanced cortical 
and subcortical thicknesses: Abnormally thin cortices grow at higher 
wavelengths with increased gyrification indices (Fig. 4, left); 
Abnormally thick cortices grow at lower wavelengths with decreased 
gyrification indices (Fig. 4, right). In agreement with the literat- 
ure 18,33 , our simulations suggest that the gyral wavelength increases 
with cortical thickness. Severely thickened cortices may even sup- 
press the formation of folds entirely (Fig. 5, right). These trends are 
consistent with the pathology of diseased human brains 42 : 
Lissencephalic brains, characterized by a thick cortex (Fig. 5a) dis- 
play a smooth surface (Figs. 5b,c); Polymicrogyric brains, character- 
ized by a thin cortex (Fig. 5a) display a large number of small folds 
(Figs. 5b,c). Future studies of developing malformations in preterm 
infants would be desirable to reveal whether an abnormal cortical 
thickness is a cause or consequence of the disease. 

Our simulated quantitative metrics for brain surface morphology, 
gyral wavelengths and gyrification indices, agree well with reported 
values in the literature. Our wavelengths range from k = 3.5 mm to 
k = 10.3 mm (Figs. 3,4), which is consistent with the clinically mea- 
sured intersulcal distance of k = 8.5 mm at the end of the seventh 
month of gestation 16 . Our gyrification indices take values up to 1.9 
(Figs. 3,4,7), which is slightly lower than characteristic gyrification 
indices in healthy human brains of 2.5 to 3.0 43 . We attribute these 
discrepancies to our simplified simulation domain and to the uni- 
formity of our folding pattern. Yet, the tendency of increased gyr- 
ification 44 in regions with reduced cortical thickness 45 in 
schizophrenic brains is consistent with our simulation. 

Our current model has a few spatial and temporal limitations. 

Spatially, our most relevant limitation is isotropic volumetric 
growth 46 : In the cortex, we could refine volumetric growth as the 
collective result of two independent mechanisms, surface growth 
and surface thickening 16 . In the subcortex, we could refine growth 
as transversely isotropic, based on regionally varying axon 
orientation maps 33 . Anisotropic growth would introduce less 
regular, more directionally oriented gyri and sulci, which would 
agree more closely with the magnetic resonance images of human 



brains (Figs. 2,5,6,7). With only minor modifications, we could also 
model individual laminar growth rates across the cortical thickness 42 
and individual regional growth rates across the cortical surface 47 . 
Temporally, our most relevant limitation is associated with 
constant growth rates: Late neuronal migration, positive growth, 
takes place in the outer cortical layers to form gyri, while neuronal 
apoptosis, negative growth, takes place in the inner layers to form 
sulci 48 ' 49 . The outer layers contain sensory and associative neurons, 
which receive many axon collaterals, while the inner layers contain 
mainly projection neurons, which send one axon each 50,51 . To 
functionally link these cellular and molecular events to cortical and 
subcortical growth 37 , we could turn different growth rates on and off 
to mimic the sequence of events during gyrogenesis, including 
neuronal proliferation, differentiation, apoptosis, dendrogenesis, 
synapsogenesis, glial proliferation, lamination, and cellular 
rearrangement 9 . 

Relevance of mechanics in human brain development. We have 
established a mechanical model to explain cortical folding in the 
developing human brain. This proof of concept might have to 
undergo a few refinements, but ultimately, it may help us to 
address basic scientific questions on different scales: On the patient 
scale, can we identify candidates at risk for neurological disorders, 
and if so, how can we quantify the level of risk? On the organ scale, is 
there a general pattern of gyrogenesis for different evolutionary 
lineages, and if so, how is it correlated to growth and cortical 
thickness? On the tissue scale, is there a biological limit for gyral 
lengths, and if so, what are its key contributors? On the cellular 
scale, what are the rate limiting factors for growth, and if any, how 
can we manipulate them? Bridging the scales, is there a correlation 
between axonal metabolism and the energy of growth, and if so, what 
is the efficiency index of the living system? Ultimately, collectively 
answering these questions may bring us closer to answer the key 
clinically question in brain development: How do we systema- 
tically manipulate brain surface morphologies through pharmacolo- 
gic or therapeutic treatment? 

Conclusion. A mechanical model for brain development can explain 
variations in gyral wavelengths and gyrification indices of healthy 
and diseased human brains. A misbalance in cortical and subcortical 
growth or thickness causes morphological abnormalities: A slower 
growing or thinner cortex generally enhances folding and increases 
the gyrification index; A faster growing or thicker cortex reduces 
folding and decreases the gyrification index. These characteristics 
are in excellent agreement with the classical pathologies of 
lissencephaly and polymicrogyria. Understanding the mechanisms 
of cortical folding in the developing human brain has direct 
implications on the diagnostics and treatment of neurological 
disorders, including severe retardation, epilepsy, schizophrenia, 
and autism. 

Methods 

Continuum model for cortical and subcortical growth. Fig. 8 illustrates our model 
for cortical and subcortical growth in the developing human brain. We model growth 
using the nonlinear field theories of mechanics supplemented by the theory of finite 
growth 52 . This results in a set of five equations, which define kinematics, constitutive 
behavior, mechanical equilibrium, growth kinematics, and growth kinetics. 

Kinematically, we decompose the deformation gradient F, the spatial gradient of 
the nonlinear deformation map, into an elastic contribution F e and a growth con- 
tribution F g , 

F = F e F g . (1) 

Similarly, we decompose the overall volume change / = / e / g = det(f) into an elastic 
volume change / e = det(F) and a growth -induced volume change J s = det(F g ). 
Constitutively, only the elastic contribution F e generates stresses 41 , 

t=[X \n{f)-ii}I+iiF tl -F\ (2) 

where X and ft are the elastic Lame constants, lis the second order unit tensor, and F et 
is the transpose of F e . These stresses enter the mechanical equilibrium equation 
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Figure 8 | Continuum model for cortical and subcortical growth. The 

cortex, the gray matter, grows morphogenetically at a constant rate G c . 
Cortical growth induces subcortical deformation, which triggers 
subcortical growth. The subcortex, the white matter, grows at a stretch- 
dependent rate as G s (/ e — 7°), where G s mimics the axon elongation rate 
and f — f activates growth only, if the elastic volume stretch J e exceeds its 
baseline value 7°. The subcortical volume growth # 3 = ( V + A V) / V mimics 
the three-dimensional effect of chronic axon elongation $ = (L + AL) /L, 
which is activated when axons are stretched beyond their physiological 
limit. 
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which we solve numerically using the finite element method. To close the set of 
governing equations, we constitutively prescribe the evolution of growth. 
Kinematically, we assume that growth is purely isotropic, 

F* = $I. (4) 

This implies that the growth volume 7 s = # 3 is equivalent to the growth multiplier ,9 to 
the power of three. Kinetically, we hypothesize that cortical growth is purely mor- 
phogenetic and that subcortical growth is stretch-induced, 



,9= 



G c 



cortical growth 



G s (/ e — /° ) subcortical growth, 



(5) 



where G c and G s are the cortical and subcortical growth rates. We assume that the 
subcortex is primarily populated by axons and correlate its growth rate with the 
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Figure 9 | Gyrificaiton index to quantify cortical folding. The gyrification 
index is the ratio between the complete cortical contour length, solid line, 
and the outer cortical contour length, dashed line. It is used to quantify the 
degree of cortical folding and can be determined semi- automatically by 
contour tracing. 



experimentally well-characterized axon elongation rate 3039 . In contrast to cortical 
growth, subcortical growth is activated only, if the elastic volume stretch f exceeds its 
baseline value 7°, i.e., when axons are stretched beyond their physiological limit 31 ' 53 . 

Finite element model for cortical folding. To simulate cortical folding, we create a 
finite element model of a two-dimensional rectangular subsection of a human brain 
slice. Fig. 2 illustrates the simulation domain, which is 2 cm long, 1 cm wide, and 
0.05 cm thick, discretized by a total of 80 X 40 X 1 = 3, 200 tri-linear Ql elements 
and 19,926 degrees of freedom. We assume a plane strain state and fix the left, bottom, 
and right boundary nodes orthogonal to the boundary, but allow them to slide freely 
along the corresponding edge. We model the cortex as Neo Hookean elastic with a 
Young's modulus E c = 9,210.87 Pa and a Poisson's ratio v c = 0.458, corresponding to 
Lame constants A c = 34, 335 Pa and fi c = 3, 159 Pa 54 . We assume that the cortex is 
three times stiffer than the subcortex, which we model as Neo Hookean elastic with 
Lame constants X s = 11, 482 Pa and fi s = 1, 053 Pa. We choose the initial cortical 
thickness to t c = 0.75 mm, the cortical growth rate to G c = 2.0, the subcortical growth 
rate to G s = 0.003, and the physiological limit to 7° = l 53 . We solve the resulting 
nonlinear finite element equations using a standard Newton-Raphson method 36 . 

Gyrificaiton index to quantify cortical folding. Fig. 9 illustrates a common clinical 
metric to quantify the degree of cortical folding, the gyrification index. The 
gyrification index is the ratio between the complete cortical contour length and the 
outer cortical contour length. In humans, it increases during ontogenesis from the 
smooth initial brain surface with values close to 1.0 towards the heavily convoluted 
adult brain surface with values between 2.5 and 3.0 43 . Lissencephalic brains have 
smaller gyrification indices 2 ; Polymicrogyric brains have larger gyrification indices 3 . 
To quantify the degree of cortical folding for varying cortical growth rates and varying 
cortical thicknesses, we determine the gyrification index semi- automatically by 
contour tracing. 
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